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Abstract 

D , Within the set of generalized Skyrme models, we identify a submodel which 

i-Q ■ has both infinitely many symmetries and a Bogomolny bound which is saturated 

by infinitely many exact soliton solutions. Concretely, the submodel consists of 

*^i , the square of the baryon current and a potential term only. Further, already on the 

^. ■ classical level, this BPS Skyrme model reproduces some features of the liquid drop 

£/") ' model of nuclei. Here, we review the properties of the model and we discuss the 

CO . semiclassical quantization of the simplest Skyrmion (the nucleon). 
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1 Introduction 

The derivation of the correct low energy hadron dynamics from the underlying funda- 
mental theory (QCD) is undoubtedly one of the most prominent challenges of strong 
interaction physics, due to the non-perturbative nature of the quark and gluon interac- 
tions in the low energy limit. In the large N c limit it is known that QCD becomes equiv- 
alent to an effective theory of mesons JT), where Baryons (hadrons as well as atomic 
nuclei) appear as solitonic excitations with an identification between the baryon num- 
ber and the topological charge (2). One of the most popular realizations of this idea is 
the Skyrme model J3| i.e., a version of a phenomenological chiral Lagrangian, where 
the primary ingredients are meson fields. The original Lagrangian has the following 
form 

L = L 2 + Li + L , (1) 

where to the sigma model part 

f 2 
L 2 = -i± Tr (rfdnU U^d^U) (2) 

a quartic term, referred as Skyrme term, has to be added, to circumvent the standard 
Derrick argument for the non-existence of static solutions 

U = -^-j Tr ([tfd^U, tfd v U} 2 ). (3) 

The last term is a potential 

Lo = -fi 2 V(U,U*), (4) 

which explicitly breaks the the chiral symmetry. Its particular form is usually adjusted 
to a concrete physical situation. Finally, static properties of baryons as well as nuclei 
are derived with the help of the semiclassical quantization of the solitonic zero modes 
0,0,0. 

Within this framework, it has been established that topological solitons (Skyrmions) 
possess qualitative properties of baryons J4), 0. Quantitatively, for the original Skyrme 
model nucleon properties are reproduced within a typical precision of about 30%. For 
higher nuclei, one of the principal problems of the Skyrme model is that the binding 
energies of the nuclei result too large already at the classical level. This problem may 
be traced back to the fact that the solitons in the Skyrme model are not BPS, i.e., they 
do not saturate the Bogomolny bound of the model. Other related problems are too 
strong internuclear forces at intermediate distances and the cristal like behaviour of 
large nuclei in the standard Skyrme model. The situation seems to be better for the ex- 
citation spectra of nuclei. Indeed, by identifying the symmetries of higher Skyrmions 
and by the use of the semiclassical quantization method the (iso-) rotational excitations 
of nuclei have been described rather successfully in J6). 

In any case, if interpreted as an effective field theory, the Skyrme model should be 
viewed as a derivative expansion, where higher derivative terms have been neglected. 
Further, baryons, i.e., extended (solitonic) solutions, have regions where derivatives 
are not small, so it may not be justified to omit such terms, and one should take them 
into account. Indeed, many generalized Skyrme models, i.e., extensions of the original 
Skyrme Lagrangian to models containing higher derivative terms have been investi- 
gated 03- ED- 
Among them, the simplest and most natural generalization of the Skyrme model is 



defined by the addition of the following sextic derivative term 

Lg = ^ t Tr {<*"""& W u^d u u u^d p u)] 2 . (5) 

Phenomenologically, this term effectively appears if one considers a massive vector 
field coupled to the chiral field via the baryon density Ifl2l . IfTTl . It turns out that such 
a modified Skyrme model gives slightly better quantitative predictions Ifl2l . Ifl3ll . IfTSII . 
Ifl6l . Besides, this sextic term is at most quadratic in time derivatives, so it allows 
for a standard hamiltonian interpretation. It is, in fact, nothing but the topological (or 
baryon) current density squared, 

L 6 = A 2 7r 4 B M B^ (6) 

where 

Tr{ c ^P a U^d v UU^d p UU^d a U). (7) 
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Within the set of generalized Skyrme models, the standard Skyrme model may be 
viewed as an approximation which is mode tractable and may be sufficient for the 
description of some physical properties of nuclei. Here we propose another approxi- 
mation which consists of the potential and the sextic term only, which we will call the 
BPS Skyrme model 

L ™ = ^ [Tr (e^U^^U U^d v U U^d p U)} 2 - i^V(U, 17+). (8) 

The subindex 06 refers to the fact that in the above Lagrangian only a potential term 
without derivatives and a term sextic in derivatives are present. This model is by con- 
struction more topological in nature than any of the standard versions of the Skyrme 
model. It shares with the standard Skyrme model the stabilization of soliton solutions 
by a higher order term in derivatives (a sextic term in our case). It is an example of 
a BPS theory with topological solitons saturating the pertinent Bogomolny bound and 
therefore, with zero binding energies at the classical level. Further, this BPS model 
possesses a huge number of symmetries which lead to its integrability. What is more 
important, among its symmetries are the volume preserving diffeomorphisms on base 
space. This allows to interpret the BPS Skyrme matter as an incompressible liquid. 
Moreover, the solitons of this theory are of compact type which results (classically) 
in a finite range of interactions (contact-like interactions). As a consequence, the BPS 
Skyrme model apparently is a good guess for certain properties of nuclei. 
We want to remark that there exists another BPS generalization of the Skyrme model 
recently proposed by Sutcliffe [27] based on dimensional reduction of the (4+1) di- 
mensional Yang-Mills (YM) theory, where the SU(2) Skyrme field is accompanied by 
an infinite tower of vector and tensor mesons. 



2 The BPS Skyrme model 
2.1 Field equations 

We use the standard parametrization of U by means of a real scalar £ and a three 
component unit vector n field (r are the Pauli matrices), 



The vector field may be related to a complex scalar u by the stereographic projection 

1 



(u + u, —i(u — u), \u\ 2 — l) 



i + M 2 

For the topological current the resulting expression is 



TT 2 (l + \u\ 2 ) 4 

and, assuming for the potential V = V(tr(U + U')) = V(£) (which we assume for 
the rest of the paper) we get for the lagrangian 

Lm = ~ (i+M*y ^ vpa ^ P u a f - v 2 v(o. (9) 

The Euler-Lagrange equations read 

,.2t// n 
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the useful formulas 



d(e av P°^u p u a ) 2 



H p u» = H^u" = 0, K p e = K p u» = 0, H^ = K^u" = 2(e aVf " x ^u p u (T ) 2 . 
2.2 Symmetries 

Apart from the standard Poincare symmetries, the model has an infinite number of tar- 
get space symmetries. The sextic term alone is the square of the pullback of the volume 
form on the target space S 3 , where this target space volume form reads explicitly 

sin r 
dV = -i——±—dtidudu (10) 

(l + \u\ z y 

and the exterior (wedge) product of the differentials is understood. Therefore, the sextic 
term alone is invariant under all target space diffeomorphisms which leave this volume 
form invariant (the volume-preserving diffeomorphisms on the target S 3 ). The poten- 
tial term in general does not respect all these symmetries, but depending on the specific 
choice, it may respect a certain subgroup. Concretely, for V = V(£), the potential is 
invariant under those volume-preserving target space diffeomorphisms which do not 
change £, 

£^£, u->u(u,u,£), (1 + \u\ 2 )~ 2 d^diidu = (1 + \u\ 2 y 2 d^udu. 

Since u spans a two-sphere in target space, these transformations form a one-parameter 
family of the groups of the area-preserving diffeomorphisms on the corresponding tar- 
get space S 2 (one-parameter family because the transformations may still depend on 



£, although they act nontrivially only on u, u). Both the Poincare transformations and 
this family of area-preserving target space diffeomorphisms are symmetries of the full 
action, so they are Noether symmetries with the corresponding conserved currents. The 
latter symmetries may, in fact, be expressed in terms of the generalized integrability, as 
we briefly discuss in the next subsection. 

The energy functional for static fields has an additional group of infinitely many sym- 
metry transformations, as we want to discuss now. These symmetries are not symme- 
tries of the full action, so they are not of the Noether type, but nevertheless they are very 
interesting from a physical point of view, as we will see below. The energy functional 
for static fields reads 

E = j d * x { { iX?Y {emnllUunTll)2 + ^ V ) (U) 

and we observe that both d 3 x and € mnl i^ m u n ui are invariant under coordinate transfor- 
mations of the base space coordinates Xj which leave the volume form d 3 x invariant. 
So this energy functional has the volume-preserving diffeomorphisms on base space as 
symmetries. These symmetries are precisely the symmetries of an incompressible ideal 
fluid, which makes them especially interesting in the context of applications to nuclear 
matter. Indeed, the resulting field theory is able to reproduce some basic features of the 
liquid droplet model of nuclei, see e.g. |29l , 1371 and the discussion below. 

2.3 Integrability 

The BPS Skyrme model is integrable in the sense that there are infinitely many con- 
served charges. Indeed, it belongs to a family of models integrable in the sense of the 
generalized integrability [38], l39l . To show that we introduce 

r* 1 = — 

(l + M 2 ) 2 ' 

The currents are 

where G(u, u, £) is an arbitrary real function of its arguments. Then, 
0" J„ = Guuu^ + Gu U u^ + G u d^ - Guuu^ - G uu u^ - G u d^ 

+G a ^^ i - Gu^fj,^ = 0, 

where we used that u^K,^ = ^/C^ = 0, u^JC^ = u^JC 11 , which follow from the pre- 
vious identities. Finally using the field equations for the complex field, one, indeed, 
finds an infinite number of conserved currents. These currents are a higher dimensional 
generalization of those constructed for the pure baby Skyrme model BD1 and are gener- 
ated by the relevant subgroup of the volume-preserving diffeomorphisms on the target 
space as discussed in the previous subsection, see also RTI . 

2.4 Bogomolny bound 

In the BPS Skyrme model, there exists a Bogomolny bound for static finite energy 
solutions []. The most concise version of the proof has been given by M. Speight [], 



and here we use this version of the proof. In a first step, the energy functional has to be 
rewritten as a complete square plus a remainder, 

d * x ( n+l^pp 6 " 1 "'^""""' ± ^ 
I ,3 2/xA sin 2 £W mnl 

I d x (i + |u| 2 ) 2 e ^ mUnUl 



> ±(2\fj,ir 2 ) 



-i f 3 ^ sin 2 £_W mnl 
tH J (1 + |u| 2 ) 2 



^bps, (12) 



and in a second step it remains to prove that the remainder £bps is topological and 
does not depend on the field configuration. For this we note that the integrand in the 
above expression is not only a three-form on the base space but, in addition, it is the 
pullback under the map U of a three-form on target space. Therefore, the base space 
integral of the base space three-form is equal to the target space integral of the target 
space three-form times the number of times the target space is covered while the base 
space is covered once (i.e., times the winding number \B\), 

d^dudu sin £ 

'^J s3 (i + M 2 ) 2 



E BPS = 2x^b\^I ^—T^ VW) 



= 2A^7r 2 < y/V > 5 3 \B\ (13) 

where < \/V >$3 is the average value of \fV on the target space S 3 , 



/- J s3 dnVv 

< VV >,<?3= ^-p — — 



J s3 dn 

and dfl is the volume form on S 3 (and J„ 3 dil = tt 2 ). The corresponding Bogomolny 
(first order) equation is 

T^iL^mtlnfil = TVW (14) 

(1 + |u| 2 ) 2 
and is satisfied by all soliton solutions which we will encounter below. 

2.5 Exact solutions 

We are interested in static topologically nontrivial solutions. Thus u must cover the 
whole complex plane (n covers at least once S 2 ) and £ <G [0, tt] . The natural (hedgehog) 
ansatz is 

e = £(r), u(0,<f>) = g{0)e* n *. (15) 

Then, the field equation for u reads 

1 a ( 9 2 9e \ gg 2 



(1 +g 2 ) 2 sinfl J (1 + 5 2)2 sin 2 q 



(17) 



and the solution with the right boundary condition is 

g{6)= tan-. (16) 

Observe that this solution holds for all values of n. The equation for the real scalar 
field is 

n 2 X 2 sin 2 £ f sm 2 ^ r \ , 

2r2 * (— ^g— j - M V{ = 0. 

This equation can be simplified by introducing the new variable 

_ V2>r 3 
Z ~ 3|n|A ' 

It reads 

sin 2 £ 9, (sin 2 ££ z ) ~ Vi = 0, (18) 

and may be integrated to 

I S m 4 ^ 2 = ^ (19) 

where we chose a vanishing integration constant to get finite energy solutions. We 
remark that this first integration of the field equation is just the dimensionally reduced 
Bogomolny equation. 

For a further evaluation we have to specify a potential. Here we shall restrict ourselves 
to the case of the standard Skyrme potential 

V = -Tr(l - [/) -)■ V(f) = l-cos£. (20) 

Imposing the boundary conditions for topologically non-trivial solutions we get 

^ ( 2arccos^ z € [0, |] (21) 

\ z>i 

The solution is of the compacton type, i.e., it has a finite support [43] (compact solu- 
tions of a similar type in different versions of the baby Skyrme models have been found 
in [33), H4l ). The corresponding energy is 

E = 8\/27r^A|n| / [ 1- ( -^ J 3 ) dz= ~ n\\n\ (22) 

and is linear in the topological charge \B\ = \n\. The energy density and topological 
charge density per unit volume are 

£ = 8-\/2juA(1- \n\~ir 2 ) for < f < \n\^ 

= for f>\n\i. (23) 

B = sign(n) — (1 — \n\^^f 2 )^ for 0<f<|n|3 

= for f>|n|3 (24) 

where 

(here Rq is the compacton radius for the soliton with charge one), see Figure 1 . 
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Figure 1 : Normalized energy density (left figure) and topological charge density (right 
figure) as a function of the rescaled radius f, for topological charge n=l. For \n\ > 1, 
the height of the densities remains the same, whereas their radius grows like \n\ 3 



3 Some phenomenology of nuclei 

For simplicity, we choose the standard Skyrme potential V 
section. 



1 — cos £ throughout this 



3.1 Classical aspects 

Already the classical solutions of the BPS Skyrme model describe quite well some 
static properties of nuclei |f29| . Concretely, we find 



• A linear mass - baryon number relation: E = 64"\/2/iA|i?|/15 = £" 1 i3 1 , which 
is a well-known relation for nuclei 

• No binding energies due to the BPS property of soliton solutions. Binding ener- 
gies for physical nuclei are small (< 1%), while in the standard Skyrme model 
they are significantly bigger. 

• No forces between solitons: due to the BPS and compact nature of solutions there 
are no long range forces. Physical nuclei have a very short range interaction. 

• Radii of nuclei: the compact solutions have a definite radius R = (2y / 2A|i?|//.j)^ 1 / 3 ^ 
= Rq \J\B\ which, again, is a well-known relation for physical nuclei. 

• Incompressible fluid: the static energy functional has the volume-preserving dif- 
feos (VPDs) as symmetries, which are the symmetries of an incompressible ideal 
fluid. Physical nuclei do not have this symmetry: they have definite shape, and 
deformations cost energy. But: volume-preserving deformations cost much less 
energy, as a consequence of the liquid drop model of nuclei. 

To summarize, the classical model already reproduces some features of the nuclear 
liquid drop model (mass and volume proportional to baryon charge, strictly finite size, 
VPDs). But clearly this picture is too naive (there are no pions and, therefore, no long 
range interactions and no pion cloud; there are no quantum corrections, etc.). Further, 
these classical results are probably more trustworthy for not too small nuclei, because 

i) the contribution of the pion cloud (which is absent in our model) to the size of 
the nucleus is of lesser significance for larger nuclei. We remind that in addition 



to the core of a nucleus (with a size which grows essentially with the third root 
of the baryon number) a surface term is known to exist for physical nuclei whose 
thickness is essentially independent of the baryon number. 

ii) the description of a nucleus as a liquid drop of nuclear matter is more appropriate 
for larger baryon number. 

iii) the contribution of (iso-) rotational quantum excitations to the total mass of a 
nucleus is smaller for larger nuclei, essentially because of the larger moments of 
inertia of larger nuclei. 

3.2 Quantization 

Here we perform the first steps towards a semiclassical (collective coordinates) quanti- 
zation of a soliton, as is done for the standard Skyrme model. Concretely, the standard 
procedure consists in quantizing the (iso-)rotational degrees of freedom of each soli- 
ton. In principle, one can quantize all solitons in this way and obtain information about 
the binding energies and rotational spectra of nuclei, but here we shall restrict to the 
simplest B — 1 case (the nucleon). Following the standard treatment, we introduce the 
collective coordinates of the isospin by including a time-dependent iso-rotation of the 
classical soliton configuration 

U{x) = A{t)U {x)A\t), (26) 

where A(t) = ao + ia^ji <G SU(2) and a^ + a 2 = 1. We remark that the B = 1 hedge- 
hog is invariant under a combined rotation in base and isotopic space, therefore, it is 
enough to introduce the collective coordinates of one of the two. The standard proce- 
dure consists in inserting the iso-rotation into the lagrangian, introducing the conjugate 
momenta to the SU(2) variables a n and applying the standard canonical quantization 
to the latter. This results in an energy expression which is the sum of the static soliton 
energy E and an isospin (or equivalently spin) contribution 

h 2 i 2 h 2 s 2 

H — Eq H — — — = Eq 
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where I 2 is the isospin operator (the spherical laplacian on S 3 ). Further, 1 is the 
moment of inertia 

T A7T I A ■ 4 fC 2 A/27T (X 

The soliton with baryon number one is quantized as a fermion. Concretely, the nucleon 
has spin and isospin 1/2, whereas the A resonance has spin and isospin 3/2, so we find 
for their masses 

Q/j2 icfc2 ofc2 

M N = E + —, M A = E + — =*> M A -M N =—, (27) 

like in the standard Skyrme model. These expressions may now be fitted to the physical 
masses of the nucleon (M N = 938.9 MeV) and the A resonance (M A = 1232 MeV), 
which determines the coupling constants 

Xu = 45.70McV , - = 0.2556fm 3 



radius 


BPS 


massive 


exp. 




Skyrme 


Skyrme 




compacton 


0.90 


- 


- 


electric isoscalar r e ,o 


0.64 


0.68 


0.72 


electric isovector r e ,i 


0.67 


1.04 


0.88 


magnetic isoscalar r m> o 


0.71 


0.95 


0.81 


?"e,l/ r e,0 


1.05 


1.53 


1.22 


r m ,o/r e ,o 


1.12 


1.40 


1.13 


r e ,i/r m ,o 


0.94 


1.10 


1.09 



Table 1 : Compacton radius and some charge radii and their ratios for the nucleon. The 
numbers for the massive Skyrme model are from Ref. (5). All radii are in fm. 



These may now be used to "predict" further physical quantities like, e.g. the charge 
radii of the nucleons. For the isoscalar (baryon) charge density per unit r we find 



p Q =4TTr 2 B° = sin 2 £C 

and for the isovector charge density per unit r 

P\ = jyA sin ££ r . 
The corresponding isoscalar and isovector mean square electric radii are 



< r 2 > e , = / drr 2 p = - 



< r' > e i = 



2/3 



drr p\ 



10 /A 



2/:>> 



9 \n 
Further, the isoscalar magnetic radius is defined as the ratio 



2 / drr 4 p 

< r 2 > m n= J 



/ drr 2 po 



2/3 



(28) 



(29) 



(30) 



(31) 



(32) 



The numerical values are displayed in TableQ] 

Firstly, in the BPS Skyrme model all radii are bound by the compacton radius Rq = 
\/2(\/ '/i)' 1 / 3 '. Secondly, all radii in the BPS Skyrme model are significantly smaller 
than their physical values, as well as significantly smaller than the values predicted in 
the standard massive Skyrme model. This has to be expected, because we know already 
that the pion cloud is absent in the BPS Skyrme model. Thirdly, we see that the error 
due to the absence of the pion cloud in the model partly cancels in the ratios, as one 
would expect. 

4 Conclusions 

In this work we proposed a submodel within the set of generalized Skyrme models 
which we called the BPS Skyrme model. The model is quite interesting on its own, 
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because of its infinitely many symmetries and conserved charges charges, its solvability 
for any form of the potential, and because all solutions are of the BPS type. They obey 
a first order differential equation and saturate a Bogomolny bound. 
Further, we found some evidence that the model may be a good starting point for an 
effective field theory description of nuclei. For a further development of tis application 
of the model, however, the following problems have to be resolved or investigated. 

• Symmetry breaking: the infinitely many symmetries of the model are not shared 
by physical nuclei. In addition, it is not clear how to select the correct soliton 
from the infinitely many ones related by symmetries or how to quantize these 
infinitely many symmetries. Therefore, a realistic phenomenological application 
will require the breaking of these symmetries. The challenge will be to identify 
a breaking mechanism which effectively breaks the unwanted symmetries with- 
out perturbing too much the good properties of the model (like weak binding 
energies, weak internuclear forces). 

• Quantization of higher nuclei: the semiclassical quantization of higher solitons 
should be performed and applied to higher nuclei, along the lines of what was 
done for the standard Skyrme model, e.g., in (6). Recently in J49l . the authors 
used a version of the BPS Skyrme model with a different potential for this pur- 
pose. Concretely, they calculated the exact static soliton solutions plus the (iso- 
) rotational energies in the rigid rotor quantization for general baryon charge 
B = n for the spherically symmetric ansatz. Then they allowed for small con- 
tributions to the total energies from the quadratic and quartic Skyrme terms and 
fitted the resulting binding energies to the experimental binding energies of the 
most abundant isotopes of higher nuclei, assuming, as is usually done, that these 
correspond to the states with the lowest possible value of the isospin. The result- 
ing agreement between calculated and experimentally determined masses and 
binding energies is impressive, lending further support to the viability of the BPS 
Skyrme model as the leading contribution to an effective theory for the proper- 
ties of nuclear matter. Their specific choice of the potential was motivated by 
the two requirements to avoid the compact nature of the solitons, and to have the 
standard pion tail in the full model with the quadratic Skyrme term included. It 
is this point where we slightly disagree, because we think that these two require- 
ments are mutually incompatible (this issue is, however, completely irrelevant 
for the remaining results of that paper). Their potential has, in fact, two vacua 
and behaves quadratically about one vacuum, but sextic around the other. Fur- 
ther, their solitons approach the sextic vacuum at large distances. This implies 
that the pure BPS model solutions are not compact, but in the full model there 
is no pion tail (the pion tail would be induced by the quadratic vacuum which 
is, however, assumed at the center of the soliton, and not in the asymptotic large 
distance region). In any case, further investigations of the semiclassical quanti- 
zation with the inclusion of higher excitations of nuclei and electrostatic effects, 
among others, would be very interesting. 

• Motivation from QCD: it would be very interesting to see whether the rather 
good phenomenological properties of the model can be justified in a more rigor- 
ous manner from the fundamental theory of strong interactions, i.e., QDC. In this 
context we want to emphasize that the two term of the BPS Skyrme model are 
rather specific. The sextic term is essentially topological in nature and, therefore, 
naturally related to collective excitations of degrees of freedom of the underlying 
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microscopic theory. The potential, on the other hand, provides the chiral sym- 
metry breaking and might therefore be related to collective degrees of freedom 
of the quarks, like the quark condensate. These issues certainly require further 
investigation. 

In any case, we think that we have found and solved an important submodel in the space 
of generalized Skyrme models, which is singled out both by its capacity to reproduce 
qualitative properties of the liquid drop model of nuclei and by its unique mathematical 
structure. The model directly relates the nuclear mass to the topological charge, and 
it naturally provides both a finite size for the nuclei and the liquid drop behaviour, 
which probably is not easy to get from an effective field theory. Finally, our exact 
solutions might provide a calibration for the demanding numerical computations in 
physical applications of more general Skyrme models. 
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